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STABILITY  REGIMES  IN  A  ROTATING 
QUADRUPOLE  FOCUSING  ACCELERATOR 

I.  INTRODUCTION 

A  number  of  recent  high  current  accelerator  configurations  utilize  strong  focusing 
fields.  These  fields,  consisting  of  a  stellarator  field  (or  rotating  quadrupole  field)  and  an 
axial  guide  field,  increase  considerably  the  energy  mismatch  tolerance  of  the  device  and 
provide  confining  forces  against  the  beam  space  charge  forces.1,2  Two  such  devices  are 
the  modified  betatron  accelerator3  and  the  spiral  line  induction  accelerator  (SLIA).4-9 
The  recent  addition  of  strong  focusing  to  the  modified  betatron  at  the  Naval  Research 
Laboratory  has  allowed  that  accelerator  to  successfully  accelerate  a  0.5  kA  beam  to  12 
MeV  for  cz  35,000  turns  in  a  vacuum  chamber  with  resistive  walls.7  The  SLIA,  to  be 
constructed  by  Pulse  Sciences,  Inc.,  will  utilize  the  strong  focusing  for  transport  along  the 
curved  sections  of  the  beam  line  between  the  accelerating  cavities. 

The  use  of  strong  focusing  fields  has  a  potential  difficulty  in  that  they  can  lead  to 
various  types  of  beam  instabilities.8,9  It  has  been  suggested,  for  example,  that  the  rotating 
quadrupole  field  may  act  like  a  wiggler  field  as  in  the  free  electron  laser.8’10  Of  interest  in 
this  study,  however,  is  the  three-wave  instability  in  which  the  externally  imposed  rotating 
quadrupole  field  interacts  with  the  transverse  motion  of  the  beam  centroid  to  excite  a 
transverse-electric  (TE)  waveguide  mode.9  Approximate  stability  conditions  for  such  strong 
focusing  fields  have  recently  been  found.11 

The  purpose  of  the  present  study  is  to  obtain  detailed  analytical  conditions  for  stabiliz¬ 
ing  the  three- wave  instability.  Here  we  analyze  the  dynamics  of  electromagnetic  waveguide 
modes  and  the  beam  modes  associated  w'  n  a  relativistic  electron  beam  propagating  under 
the  influence  of  external  strong  focusing  fields.  The  dispersion  relation  for  such  a  system 
is  generated  in  Section  II  below.  In  Section  III,  we  show  that  this  dispersion  relation,  in 
the  limit  of  zero  beam  current,  gives  analytical  conditions  for  the  various  stability  regimes. 
We  find  a)  two  physically  distinct  three-wave  unstable  regimes,  b)  two  three-wave  stable 
regimes  and  c)  a  regime  in  which  the  particle  orbits  themselves  are  unstable,  irrespec¬ 
tive  of  the  electromagnetic  waves.  The  orbit  unstable  regime  is  a  recovery  of  an  earlier 
result.1,2  In  addition,  we  obtain  algebraic  expressions  for  the  growth  rates  in  each  of  the 
three-wave  unstable  regimes.  This  is  done  with  the  simplified  dispersion  relation  where 
the  rigbUhand  (RH)  and  left-hand  (LH)  circularly  polarized  waves  are  decoupled  in  the 
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low  current  regime.  These  are  found  to  be  in  general  agreement  with  the  approximate 
expressions  given  in  Ref.  12.  Numerical  solutions  of  the  dispersion  relation  are  presented 
in  section  IV.  The  analytic  expressions  of  the  growth  rates  and  stability  boundaries  are  in 
good  agreement  with  the  numerical  results  of  the  dispersion  relation. 

II.  MODEL 

In  our  model  the  external  fields  consist  of  a  periodic  rotating  quadrupole  field  and  a 
longitudinal  magnetic  field  (see  Fig.  1).  The  electron  beam  radius  is  assumed  to  be  small 
compared  to  the  waveguide  radius.  In  the  equilibrium  position,  the  beam  is  centered  along 
the  axis  of  a  circular  waveguide  and  the  beam  electrons  are  assumed  to  be  monoener- 
getic  with  zero  transverse  velocity.  Perturbing  electromagnetic  waveguide  fields  cause  the 
beam  centroid  to  develop  a  transverse  velocity  and  become  displaced  off  the  z-axis.  This 
displacement  of  the  beam  centroid  amounts  to  a  transverse  macroscopic  current  which,  un¬ 
der  certain  conditions,  further  excites  the  electromagnetic  field.  The  displacement  of  the 
beam  centroid  also  induces  image  electric  and  magnetic  fields  on  the  wall  of  the  waveguide. 
The  beam  centroid  motion,  in  our  model,  is  governed  by  the  following  fields:  i)  rotating 
quadrupole  field,  ii)  longitudinal  magnetic  field,  iii)  electromagnetic  fields,  and  iv)  induced 
image  fields. 


Wave  equation 

Before  developing  the  orbit  equations  for  the  beam  centroid  we  first  derive  the  wave 
equation  for  the  electromagnetic  fields.  The  electromagnetic  fields  are  represented  by  a 
vector  potential  given  by, 

A  =  A+(x,  y,  z)e~lwt  +  A_(z,  y,  z)e~,u,t  +  c.c.,  ( 1 ) 


where  A+  and  A_  are  complex  amplitudes  associated  with  the  right-hand  (RH)  and  left- 
hand  (LH)  circularly  polarized  waves  in  a  cylindrical  waveguide.  In  Eq.  ( 1 ),  u>  is  the  radian 
frequency  and  c.c.  denotes  the  complex  conjugate. 

The  wave  equation  for  A  is  given  by 


/_2  d2  1  d2  \  A  4tt 

V1  +  dz2  c2dt2  j  A  ~  c  J’ 


(2) 
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where  Vj_  is  the  transverse  Laplacian,  and  J  is  the  macroscopic  transverse  current  asso¬ 
ciated  with  the  beam  centroid.  Substituting  (1)  into  (2)  and  operating  on  the  result  with 
(<±>/2t)  J'**'  u  dt  exp(iu/t),  in  order  to  select  the  correct  frequency  dependence,  yields 


_  ,  .  .  4 7T  f2*1"  dt  •  . 

I(A++A_)  = - /  -  re  J, 

c  J0  2tt j  uj 


(3) 


where  L  =  -I-  d2 /dz*  -fa ;2/c2.  To  obtain  the  RH  and  LH  polarized  components  of  the 

current,  we  equate  the  e+  and  e_  components  of  both  sides  of  (3),  and  find  the  following 
wave  equation, 


4ir  r2*/"  dt 

L(A+  +  A_)  •  (2eT)  - - /  —  (Jcx  T  iJcV)eluJ\ 

c  y o  2irjui 

where  e±  =  (ex  ±  iey)/ 2  and  ex  and  ey  denote  unit  vectors  in  the  x-  and  y-directions. 
The  current  density  associated  with  the  beam  centroid  motion  is 


(4) 


J{x,y,z,t)  =  -  |e|A(, 


/>( 


dxc(z0,t) .  _  dyc(z„,t)„ 

~ — e*  +  ~dT~^ 


(5) 


5(i  -  xc(za,t))6(y  -  yc(z0,t))S(z  -  ze(z0,t)), 


where  A<,  =  nbirrl  is  the  number  of  electrons  per  unit  length,  nb  is  the  beam  density,  rb 
is  the  beam  radius,  xe  and  yc  denote  the  transverse  coordinates  of  the  beam  centroid. 
zc  denotes  the  axial  position  of  a  cross  sectional  slice  of  the  beam  and  z„  is  the  initial 
position  of  the  slice:  zc(z0,t  =  0)  =  za.  In  the  small-signal  or  linear  regime  xc  and  yc  are 
proportional  to  the  electromagnetic  fields  and  (5)  can  be  written  as 


Jc(x,y,z,t)  ~  -  |e|Aj,<5(x)£(y) 


/X^r 


,  9yc(zo,0- 


ex  -f 


dt 


ey  )  S(: 


:C  (ZOl  t)), 


(6) 


where  v0  is  the  axial  beam  velocity.  Substituting  (6)  into  (4),  we  obtain 

L(  A+  +  A_)  •  (2eT)  =-u6(x)6(y)  f  dza 

c  JQ  2tt/ u>  J_<x 


dxc(z0,t)  .dyc{z0,t) 

T"  * 


dt 


dt 


j  6(z  -  z0  -  v0t)el 


~) 
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where  u  =  (|e|/m0c2)/&//?0  ~  Ib[kA\/\7f30  is  Budker’s  parameter,  It,  is  the  beam  current 
and  (30  =  v0/c.  Upon  carrying  out  the  time  integration  in  (7)  we  obtain, 


L(A+  +  A_)  •  (2eT)  =  —  u6(x)6{y)el{u'/v'’)z 

c 


L 


2wv„/u  i  o 

■+■  lyc,ie  j 

o  Zttv0/(jJ  ai 


(8) 


where  xc(z0,z)  —  xc (zQ,t  —  (z  z0'j/v0'j  and  ye(za,z)  —  yc{z0,t  —  (z  z0)/u„).  Since 

our  model  is  spatially  periodic  with  period  u>/2ir,  the  limits  on  the  z„  integral  have  been 
changed  to  0  to  2ir/u. 


Beam  Centroid  Orbit 

The  configuration  of  the  rotating  quadrupole  and  the  axial  magnetic  fields  are  shown 
in  Fig.  1.  The  total  external  magnetic  field,  Bej.t,  consists  of  the  rotating  quadrupole  field 
and  the  uniform  longitudinal  field  and  is  given  by  Be*t  =  (Bqx,  Bqy ,  Bzo)  where 

Bqx  =  -Bqkq(x  sin  kqz  -  y  cos  kqz), 

Bqy  =  Bqkq(x  cos  kqz  +  y  sin  kqz),  (0 a  -  c) 

Bzo  =  B0. 

In  Eqs.  (9a-c),  Ba  is  the  axial  magnetic  field,  Bq  is  the  magnetic  field  of  the  quadrupole, 
kq  =  2ir/\q  and  Xq  is  the  period  of  the  quadrupole  field.  The  representation  for  the 
quadrupole  field  in  Eqs.  (9a, b)  is  valid  near  the  z-axis,  i.e.,  (x2  +  y2)1^2  <<  A7/2t. 

In  addition  to  the  electromagnetic  fields  given  by  (1)  and  external  fields  given  by  (9), 
we  have  induced  fields  due  to  the  displaced  beam.  These  fields  are  produced  by  the  image 
charges  and  currents  on  the  waveguide  wall  when  the  beam  is  displaced  off  the  z-axis.  For 
a  circular,  perfectly  conducting  waveguide  the  induced  electric  and  magnetic  fields  near 
the  z-axis  are 


—  2  .  ,  2  "b  yc® y)i 

1*1  rg 

(10  a) 

m  c2  u 

®»n<i  —  2  .  2/io(yCe*  ■ce®y)> 

\e\  rg 

(10  b) 
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where  r3  is  the  waveguide  radius  and  we  have  assumed  (  x2c  +  yc2)1/2  <<  rg.  The  motion 
of  the  beam’s  centroid  under  the  influence  of  the  fields  in  Eqs.  (1),  (9)  and  (10)  is  in  the 
linear  approximation  governed  by 

~  v0kq£lq{xe  cos  kqz  -f  ye  sin  kqz)  -  v20k]xc 


70m0c 


( d  d 


Ax  i 


:  11a) 


4-  v0kqQq(yc  cos  kqz  -  xc  sin kqz)  -  v20k]yc 


(l+t,4)A»’ 


(116) 

70m0c  *  —  —  1  ' 

where  f 10  =  \e\B0/y0m0c  is  the  relativistic  cyclotron  frequency  associated  with  the  axial 
field,  flq  =  \e\Bq/y0m0c  \s  the  relativistic  cyclotron  frequency  associated  with  the  rotating 
quadrupole  field,  ya  =  (1  — is  the  relativistic  mass  factor,  k,  =  (2is/(/32y2r2))1/2 . 
z  =  z0-k  vat ,  and  Ax,  Ay  are  the  x  and  y  components  of  the  vector  potential  given  in  ( 1 ). 
For  the  purpose  here,  we  will  approximate  the  exact  expressions  of  (1)  by  the  fields  on 
axis, 

I  > 

■A±(x  =  0,y  =  0,  z)  -  a±elk±ze±, 


where  k±  are  the  axial  wave  numbers. 

The  orbit  equations  for  the  beam’s  centroid  can  be  written  in  a  more  convenient  form. 
Setting  £  =  xc  4-  iyc,  Eqs.  (lla,b)  become 


where  Ka  =  f20/v0,  Kq  =  ttq/v0,  and 

F  =  [-i(fc;  -u ,/v0)ale-'{k+-“/v~)ze-luz"/v-' 

Ho*1o  L 

+i(k.  -  (jj/v0)a-e'(k~  ~UJ/V°)z eiwz” /°" 

It  is  convenient  at  this  point  to  introduce  a  transformation  from  the  quantity  £  to  £ 

£  _  £'e+»fcW2_ 


12) 


14) 
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Substituting  (14)  into  (13)  transforms  the  beam’s  controid  equation  into 


r  a2 

dz2 


—  iKi 


d_ 

dz 


i  -  k-c 


—  G, 


(151 


where  K\  =  Ka  -  kq,  Af  =  (K0  —  kq/2)kq/2  —  fc2,  K$  =  Kqkq  and  G  =  Fex p(  -ikqz;2). 
Equation(15)  can  be  written  in  the  form 


jjT  +  (*  Kl  +  K\)^  +  (Ki-Ki) 


(  =  H, 


( 16a) 


where 


H  = 


0*  9! 

a?  +  'K'8*+Ki 


G  +  KiG ' 


d2  d 

— +i(A'i  -  A:,)-  4- ((A, 


kq/2)kqj2  +  A,2) 


4- A,2FVfc^/2. 


Substituting  (13)  into  (16b)  yields 


( 166) 


H  = 


Polo 


(Kl  -  <fe,/2)£>;<e“*r;2  -  (Kl  4-  kj2)K\a'_t 


-iK‘_ 


./  »-• 


~T~  ~  V2 P-a-e1*'1  -  (A+  +  /fc,/2)A'2a+e,K+£]  e1"2'"'"".  (17) 

Polo 

where  K±  —  k±  —  ujjv0  ±  kq/2  and  D±  =  Ai  T-  K\K±  —  A’|.  The  particular  solution  to 

(16a)  is 

|  =  (a;a;e",ff+2  +  Pla_e~lK- z)  e~lwz''/v'' 

4-  (a_a_elK"_  2  4-  p+a+e'K'rZ)  e"~z " " ,  ( ISa) 

where 

i 

a±  ~  —  p  ^  ( k±  —  ui / v0) / R-j., 

P±  —  —  A|oe±/ D±, 

R±  =  Kl  -  (2 A'2  4-  Af)A2  +  Kt  -  A'j  ( l86  _  n 

=  (Kl  -  (d2  +  d2))  (Kl  -  (d2  -  d;)) , 
d2  =  A,2  +  A'2/2, 

d2  =  ((A2  4-  A'i/2)2  -  (A';  -  A'j  ))1/2  . 
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Subsituting  Eqs.  (18)  together  with  (14)  into  Eq.  (8),  the  right-hand  side  of  tiie  wave 
equation  becomes 


*> 

771  C" 

Airiu30 -?—8{  x)8(yy 

|e| 

\i(K±  Tfc,/2)a±d±e1,/fiT^/3':  1191 

+(A'*  T  kq/2}fS^a^etiK^k''r-)z  . 

Dispersion  Relation 

To  obtain  the  dispersion  relation  the  coupled  differential  equation  must  be  solved, 
subject  to  the  boundary  conditions  on  the  waveguide.  To  this  end,  we  assume  that  the 
electron  beam  propagates  within  a  perfectly  conducting  cylindrical  waveguide  of  ratlins 
rg.  In  general,  the  complex  amplitude  for  the  right-  and  left-hand  circularly  polarized  TE 
waves  are  written  as, 

A±(r,0,z)e-‘“‘  =  ^ b±nm  [/n-l  (Mnmf)e±  + 

n,m 

exp  [i{k±nmz  ±  (n  -  1)0  —  ud)] ,  (20) 

where  Jn  is  the  nUi  order  Bessel  function,  b±nm  are  complex  constants.  are  real 

constants  determined  by  the  boundary  conditions  and  n  =  1, 2, 3. ...  and  m  =  1,2.3....  are 
waveguide  mode  indicies. 

The  boundary  condition  is  such  that  the  tangential  component  of  the  electric  field 
vanishes  on  the  waveguide  surface,  r  =  rg.  At  r  =  rg  we  have  E  •  e#  =  0,  where  E  is 
the  total  electric  field  and  =  —  iexp(  —  i’0)e+  -f  iexp(  —  i#)e_  is  the  unit  vector  in  the 
azimuthal  direction.  Applying  the  boundary  condition  at  r  =  rg,  we  find  the  condition 
J'n(l*nmrg)  =  30  that  iinmrg  equals  the  mth  positive  zero  of  J'n. 

Since  A  is  driven  by  an  effective  transverse  line  current,  we  expect  TE i„,  modes  will 
be  excited.  Substituting  the  TE\m  mode  representation  into  the  left-hand  side  of  Eq.  (4). 
we  obtain  two  sets  of  coupled  equations. 


47T  f 

c  Jo 


2ir/d/ 


dt 


2ir  fu) 


(Jcz  T  iJcy)e  = 
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r.  e‘*±lm*  ^  ^2  ~  ^±Im  ~  Mlm^  4(Mlmr)^±lm 

m 

V  e'*T,”’:  (^r  ~  ^lm  -  Mlrn) 


_  JiL  /2,r/w 

C  rn0c2  70 


dt 


2n/ui 


-elut(JxTiJy) 


r=0 


Operating  on  Eq.  (21)  by  both  fQ  *  dO  J0{^ilmr)rdr  and 
f**  d9e*t20  Jg1  Ji(n\mr)rdT  and  solving  for  b±  i,m,  we  find  that 


!  21 ) 


u/"  .  -I  ■> 

—  -  *ilm  -  Mir 
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'±1  m  ^ 


2  1  |  e  | 

C  Am  m0c- 


2rr 


e 


ru;t 


[•L>(Mlfn**)(*^i  ^  ldy  )  |  r_  q  +  d4  ( M  l  m  7“ )  (  -^*  —  id  y 


22) 


where  /lm  =  [j£(/iimr)  4-  J22(/iimr)]  rdr  =  (Amrl  -  l)n1^ld2{filmrg). 

The  dispersion  relation  greatly  simplifies  when  only  one  waveguide  mode,  say  the 
TEu,  takes  part  in  the  interaction,  i.e.,  n  =  1,  m  =  1,  a±  =  Lu,  k±  =  k~n. 
Substituting  (19)  into  the  right-hand-side  of  (22),  we  obtain 


—  k 


±1  m 


Am 


R± 


b±n*lk±: 


=  -k\Kl  j  ~  —  T  fc,  j  b^ue 


t<  k  i 


( 23 1 


where  fc2  =  2u/^0In  =  2u^\ll{‘i0{Awr)  ~  1  )*^f (  Miirs))-  The  RH  and  LH  circularly  po¬ 
larized  waves  are  simultaneously  excited  only  when  k±  —  fc_  -  kg.  Eliminating  b±n,  the 
dispersion  relation  coupling  the  RH  and  LH  circularly  polarized  waves  becomes 


r  „  _  k;(k  -  <jj!v0  4-  kg)2  D- 
u;2/c2  -  (k  4 -kg)2  -  Mix . 


R  - 


kl(k  -  u)tv0)2D+ 
u>-/c*  -  k-  -  pi; 


=  k*4A\4 


(k  -  uj/v0  +  kq 


(fc  -  uj/i’o)2 
u;2/c2  -  k2  -  Mu  /  V^2/c2  -  (fc  4-  fc,)2  -  Mu 


(24) 
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where  k  =  k+,  k~  —  k  4-  kq,  K  =  K±  =  k  —  uj/v0  4-  kq/ 2,  D±  =  A'3  T  A’i  K  -  K\  and 
A  =  R±  —  D+D _  —  A*4  =  A-4  —  (2A|  4-  A"3) A’3  4- A*4  —  A’4.  Equation  (24)  can  be  put  into 
the  form 


feh(«  -  +  fe,)2 

u;2/c2  -  (*!  +  fc,)2  -  Mil 


fcg(A:  -  u)/r0)2 

u;2/c2  -  k2  -  n2n 


and/or 

A  =  0. 


(256) 


Equation  (25a)  agrees  with  the  dispersion  relation  in  Ref.  9  with  the  vertical  field  set 
to  zero.  Equation(25b)  is  the  dispersion  relation  of  the  hybrid  cyclotron  and  quadrupole 
modes  with  image  fields  in  a  waveguide  in  the  absence  of  electromagnetic  fields. 

The  dispersion  relation  possesses  the  following  symmetry: 


{k,kq,B0)  <!=>  (Ic  -  kq,  -B0). 

Utilizing  this  symmetry  condition,  the  discussions  in  Section  III  can  also  be  applied  to 
rotating  quadrupoles  of  the  opposite  heiicity.  The  polarization  of  the  waveguide  mode 
associated  withe  the  wave  number  k  that  the  three- wave  is  unstable  is  determined  by  the 
polarization  of  the  waveguide  mode  that  intersects  the  unstable  beam  modes. 

The  modes  taking  part  in  the  interaction  can  be  conveniently  classified  by  setting  the 
beam  current  equal  to  zero.  The  dispersion  relation,  for  zero  beam  current,  reduces  to 

[k2  -  u >2/c2  +  Mu]  [(*  +  M2  -  “Vc2  +  Mn] 

[(k  -w/v0  +  kq/2)2  -{d\  +<^)]  [( fe  —  u; / u„  +  kq/2)2  -  [d\  -  dl)}  =0,  (26) 

where 

d\  =  l-(K0(K0-kq)+k2/2), 
d2  =  l-(K;(K0  -  k,)2  +  4K;k;)'/2 . 

The  first  and  second  bracketed  terms  on  the  left-hand  side  of  Eq.  (26)  represents  the 
LH  and  RH  circularly  polarized  transverse  electric  waveguide  modes,  while  the  third  and 
fourth  terms,  A,  are  hybrid  cyclotron  and  quadrupole  modes. 
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A  diagram  of  the  dispersion  relation  in  the  zero  beam  current  limit  is  shown  in  Fig. 
2.  The  parameters  for  the  plot  are  guide  magnetic  field  B0  =  —  1  JcG,  kq  =  0.5  cm-1, 
quadrupole  gradient  Bqkq  =  200  G/cm,  7„  =  5  and  drift-tube  radius  rg  =  3  cm.  The 
curves  sure  associated  with  the  RH  and  LH  polarized  waveguide  modes  respectively.  The 
cut-off  frequency  of  the  waveguide  TEn  mode  is  ,unc.  The  two  straight  solid  ( — )  lines 
correspond  to  the  beam  lines  associated  with  the  d\  —  d\  expressions,  and  the  two  dashed 
(-  -  -)  lines  correspond  to  the  beam  lines  associated  with  the  d\  4-  d\  expression. 
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III.  STABILITY  REGIMES  AND  ANALYTICAL  EXPRESSIONS  FOR  GROWTH 
RATES 

The  dispersion  relation,  Eq.  (25a),  contains  i)  a  region  of  orbital  instability  (in  the 
absence  of  the  electromagnetic  waves),  ii)  regions  of  three- wave  instability,  and  iii)  regions 
of  stablility.  In  this  section,  we  will  obtain  the  conditions  delineating  the  various  regimes 
and  find  analytical  expressions  for  the  maximum  growth  rate  in  each  of  the  three- wave 
unstable  regions.  The  stability  diagrams  are  obtained  in  (kq,K0)  space  for  given  values  of 
7 o,  rg  and  Bqkq  in  the  limit  of  zero  beam  current. 

Orbit  Unstable  Regime 

The  expression  R  =  0  is  the  dispersion  relation  for  the  particle  dynamics  in  the  pres¬ 
ence  of  stellarator  windings  with  an  axial  magnetic  field.  This  expression  is  in  agreement 
with  Eq.  (10)  of  Ref.  2  in  the  limit  of  perfectly  conducting  walls.  The  electron  beam  in 
this  configuration  can  be  unstable  when  (d\  —  d\)  <  0.  The  unstable  values  of  Ka  are 

Kcrit,2  =  ~  <  K0  <  KCrit,3  =  +  2Aq,  (27) 

where  Kq  =  Uq/va,  ttq  =  |e|B,/70m0c,  K„  =  f l„/v0  and  =  \e\B0j^omoc.  Equation  (27) 
is  in  agreement  with  the  stability  condition  of  Ref.  1  in  the  limit  of  straight  cyliudrici  ’ 
geometry  and  zero  beam  current.  It  is  interesting  to  note  that  this  condition  is  also  in 
agreement  with  the  condition  for  beam  envelope  stability  in  the  limit  of  zero  space  charge. 13 

Three- Wave  Unstable  Regimes 

The  three-wave  instability  will  occur  when,  for  example,  the  RHCP  waveguide  mode 
intersects,  in  the  (w,k)  plane,  the  appropriate  beam  mode  given  by 


(fc  -  ui/v0  -(-  kq/2)2  -  {d\  -  d\ )  =  0,  (2S) 

and  {d\  —  d\)  >  0.  For  kq  >  0,  the  instability  for  the  RHCP  waveguide  mode  occurs  for 
ijJ  >  0.  For  K0  <  KCrit,2i  the  three- wave  is  unstable  (Region  I)  when  the  RHCP  waveguide 
mode  intersects  the  beam  line  u/v0  =  (k  +  kq/ 2)  -f  \J d\  —  d\ .  For  K0  >  KcrU,z,  the  three- 
wave  is  unstable  (Region  II)  when  the  RHCP  waveguide  mode  intersects  the  beam  line 
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ui/v0  =  (k  +  kq/2)  —  y  d\  —  d\.  Identical  three- wave  instability  growth  rates  occur  for  the 
LHCP  waveguide  mode,  for  u//c  with  sign  opposite  to  the  unstable  modes  associated  with 
the  RHCP  waveguide  mode. 

Three- Wave  Stable  Regime  for  K„  <  A'-,.,, 

Stability  is  achieved  when  the  waveguide  cut  off  frequency  Mnc  is  sufficiently  large  so 
that  intersection  with  either  of  the  beam  lines,  defined  by  Eq.  (28)  cannot  be  achieved. 
The  condition  in  terms  of  the  waveguide  mode  cutoff  is 

<ZMu  >  kq  +  2  (d\  -  d\)  !  ,  (29) 

where 


Based  on  Eq.  (29),  the  region  of 

kq  >  qpi i  and  Ka  <  KCrit,2  (30) 

is  alway  three- wave  unstable.  If  the  inequality  in  (29)  can  be  satisfied,  we  can  solve  for 
the  explicit  value  of  K0  for  the  three-wave  stable  regime. 

For  kq  <  qfxu  and  K0  <  Kerit, 2,  condition  for  stability  in  terms  of  the  variable  A'0  is 

(31) 


(32) 


(.331 


where 


2fK0{K0  -k ,)  4-  /2  -  4A*,Jfc,  <  0, 


/  =  q»u{kq  -  qn\x! 2). 


Defining 

C  =  ^(l+8A-J/)-2/. 

we  solve  for  Ka  with  kq  <  q/iu,  and  find  four  situations: 
i)  for  /  >  0  and  £  >  0,  the  stable  range  of  Ka  is  given  by 


K  — 

—  2 


c 


1/2 


<  K0  <  A 


crt£,2  1 


( 34fl ) 


ii)  for  /  <  0  and  (  >  0,  the  stable  values  of  Ka  are 


Ka  <  A'cr»t,i  =  smaller 


<1/2  K 

2  >-^crtl,2 


iii)  for  /  >  0  and  £  <  0,  all  values  of 


(346) 


K0  <  K 


crit,  2 


(34  c) 


are  unstable, 

iv)  for  /  <  0  and  £  <  0,  all  values  of 


A'  <  K, 


crit, 2 


(34  d) 


are  stable. 

Three- Wave  Stable  Regime  for  K„ 

The  three-wave  interaction  is  also  stable  when  the  RHCP  waveguide  mode  intersects 
only  the  top  (the  beam  line  with  larger  u;  for  the  same  k)  but  not  the  bottom  beam  mode 
assocated  with  Eq.  (28),  and  Ka  >  KCrit,s ■  This  occurs  when 

<7^ii  >  *,  ~2{d\  -dl)l/2.  (35) 


The  three- wave  interaction  is  stable  for 

fc,  <  9Mn  and  Ko  >  K crit, 3-  (36) 

For  kq  >  qnn  and  K0  >  KeT~^t,3^  the  values  of  K„  that  are  three- wave  stable  are 

k  Cl/2 

A o  >  Acrjt,4  =  H — ~ ■  (37) 

In  this  regime,  one  can  show  that  /  >  0  and  £  >  0.  In  the  limit  of  small  quadrupole 
gradient,  luge  7<,  and  Ka  >  Kcrit,3,  the  stability  condition  is  approximately11 


Ko  >  kq  Pll/fo' 


(38) 
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The  various  operating  regimes  are  illustrated  as  functions  of  kq  and  Ka  in  Fig.  3,  for 
j0  =  5,  rg  =  3  cm  and  quadrupole  gradient  Bqkq  =  200  G/cm.  This  plot  assumes  that 
the  quadrupole  gradient  is  a  constant,  and  (kq,  K0)  are  allowed  to  vary.  The  horizontal 
separation  of  the  boundaries  for  the  orbital  unstable  region  is  4Kq.  Since  the  stability 
boundaries  are  obtained  in  the  limit  of  zero  beam  current,  the  area  of  the  actual  stable 
regions  will  shrink  slightly  as  the  value  of  the  current  is  increased. 

Analytical  Expressions  for  the  Growth  Rates 

In  the  three-wave  unstable  regimes,  we  can  obtain  analytical  expressions  for  the  peak 
growth  rates.  The  dispersion  relation  (25a)  for  the  coupled  RH  and  LH  polarized  wave 
can  be  rewritten  as 

(a/2  —  ar2)((a >  —  a ;2)2  —  Awf)((w  —  u>2)2  —  Au>|)(  ui2  —  u>\)  —  a,  (39) 


where 

9  =  kbvlc2  [(w  ~  v0k)2{uj2  -  wl)D+  +  (a>  -  va(k  +  k,))2(u;2  -u>2)D_]  , 

“>i  —  y/k2  +  MiiC,  u/2  =  v„(k  4-  kq/ 2),  u>3  =  y/(k  +  kq)2  +  jijjC,  Awj  =  vl(d]  +  d;),  and 
Ao/|  —  v2(dj  —  d2).  Based  on  the  numerical  results  of  the  full  dispersion  relation  in  Eq. 
(25a),  instability  occurs  at  the  intersection  of  the  RHCP  waveguide  mode  and  one  of  the 
two  modes  of  Eq.  (28)  for  fc,  >  0  and  w  >  0.  Defining  u>  =  u>i  4-  6w  and  <r  =  ?|w=wi,  the 
dispersion  relation  reduces  to 

(u>2  -  W2)((w  -  u>2)2  -  Au/2)((u/  -  u/2)2  -  Au/j)  =  <r,  (40) 

where  <r  =  k\v20c 2  (wt  -  v0k)2D+  D+  =  (Au;2/c)2  ± h\(  Aw2/c)  -  A',2 

and  the  top  and  bottom  signs  in  D+  refer  to  the  three-wave  unstable  regions  I  and  II 
respectively.  Equation  (40)  is  the  dispersion  relation  when  the  RHCP  waveguide  mode  is 
uncoupled  to  the  LHCP  waveguide  mode. 

The  instability  region  I  in  Fig.  3  is  the  result  of  the  waveguide  mode  intersecting  the 
upper  mode  (the  beam  line  with  larger  u  for  the  same  k)  given  in  Eq.  (28).  We  will  write 
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uj  =  ui\  -\-  8w  and  u>i  =  w2  4-  Au>2.  We  will  assume  Ao/i  —  Au>2  >>  8u.  The  dispersion 
relation  is  given  approximately  by 


<rjc3 _ 1 _ 

2u>\  (Au>j  —  Aw|) 


(41) 


For  simplicity,  we  will  assume  the  temporal  growth  rate  is  much  smaller  than  the  separation 
of  the  beam  modes,  i.e.,  6u  <<  2Aw2.  With  this  assumption,  we  obtain  a  simple  expression 
for  the  dispersion  relation, 


(Su ^/c2  1  (42) 

\  c  /  4u;1Au»2  (Au>j  —  Au;|) 

For  values  of  K0  in  the  unstable  region  I,  i.e.,  K0  <  Kcritt2,  the  values  of  cr  are  positive, 
and  Eq.  (42)  gives  the  temporal  growth  rate. 

When  the  RHCP  waveguide  mode  intersects  both  modes  given  by  Eq.  (28)  with 
kq  >  0,  u;  >  0  and  K0  >  3,  the  instability  occurs  only  at  the  intersection  of  the 

RHCP  wave  mode  and  the  lower  beam  mode.  This  gives  the  instability  regime  II  in  Fig. 
3.  We  will  define  u>i  =  w2  —  Au»2.  The  dispersion  relation  becomes 


/  6ui  Au;2\/&*;\2_  <r/c 3  1 

\  c  c  )\c  )  2u/j  (Au/f  -  A  u/|)' 

Here  again,  we  assume  2Au>2  >>  8u>  and  the  dispersion  relation  reduces  to 

/  6w\2  ^  <r/c2  1 

\  c  /  4wi  Aw2  (Awf  —  Aw|) 


(43) 


(44) 


For  Kerit,  3  <  Ka  <  kq ,  the  quantity  cr  is  negative,  and  Eq.  (44)  gives  the  temporal  growth 
rate  in  region  II  of  Fig.  3.  The  analytical  expressions  (42)  and  (44)  show  that  the  temporal 
growth  rate  T  =  /m(£u>)  scales  approximately  as  the  square  root  of  the  beam  current,  i.e., 

r/c  a  ri/2. 
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rv.  NUMERICAL  RESULTS 

The  full  dispersion  relation,  Eq.  (25a)  is  solved  numerically  to:  i)  obtain  the  growth 
rates  and  group  velocities,  ii)  verify  the  various  operating  regimes  and  the  analytical  expres¬ 
sions  for  the  temporal  growth  rates  and  iii)  show  the  scaling  of  the  growth  rate  with  respect 
to  the  various  parameters.  The  numerical  studies  center  around  the  parameters  shown  in 
Table  I.  With  the  quadrupole  wave  number  kq  chosen  to  be  0.5  cm-1  (A,  =  12.57  cm), 
we  may  demonstrate  each  of  the  different  operating  regimes  by  varying  the  axial  magnetic 
field,  B„ ,  except  in  one  of  the  stable  regimes,  where  we  take  kq  =0.1  cm'1. 

Figure  3  is  a  plot  of  the  various  operating  regimes  in  the  parameter  space  of  A:,,  as  a 
function  of  Ka  for  7<>  =  5,  Bqkq  =  200  G/cm  and  rg  =  3  cm  as  in  Table  I. 

Numerical  Results  from  Full  Dispersion  Relation 

The  dispersion  diagram  with  current  /(,  =  1  kA  for  the  five  different  regimes  are  shown 
in  Fig.  4: 

a)  Ba  =  — 1.0  kG  (K0  =  —0.24  cm'1)  and  kq  =  0.5  cm'1  in  the  three- wave  unstable 

region  I, 

b)  B„  —  2.15  kG  ( Ka  =  0.26  cm'1)  and  kq  =  0.5  cm'1  in  the  orbit  unstable  regime, 

c)  Ba  =  3.5  kG  (K„  =  0.42  cm-1)  and  kq  =  0.5  cm'1  in  the  three- wave  unstable  region 

II, 

d)  B„  =  5.0  kG  (Ka  =  0.60  cm'1)  and  kq  =  0.5  cm'1  in  the  three-wave  stable  regime, 

e)  B0  =  —5.0  kG  ( Ka  =  —0.60  cm'1)  and  kq  =  0.1  cm'1  in  the  three- wave  stable 

regime. 

Plots  of  the  temporal  growth  rate  as  a  function  of  wave  number  k  are  given  in  Figs.  5-7 
for  each  of  the  different  regimes  with  current  A  =  1  kA.  Figure  5  plots  the  temporal  growth 
rates  in  region  I,  for  Ba  =  —1.0  kG  ( K0  =  —0.12  cm'1),  B0  =  0,  B0  =  1.0  kG  (K„  = 
0.12  cm'1),  and  B0  =  1.3  kG  (K„  =  0.156  cm'1).  Only  the  growth  rates  associated 
with  JRe(w)  >  0  are  plotted;  these  are  associated  with  RHCP  waves.  The  growth  rates 
associated  with  Re( u>)  <  0  are  identical  and  are  associated  with  LHC’P  waves.  The  growth 
rate  and  the  range  of  unstable  values  of  k  increase  as  Ka  approaches  the  orbit  unstable 
value  of  Keru,2  =  0.154  cm'1. 


16 


Figure  6  shows  temporal  growth  rates  in  the  three-wave  unstable  region  II,  for  B0  = 
2.95  kG  (Ka  =  0.35  cm'1),  B0  =  3.0  kG  ( K„  =  0.36  cm'1),  B0  =  3.25  kG  (K„  = 
0.39  cm"1)  and  B„  =  3.5  kG  ( Ka  =  0.41  cm-1).  For  some  values  of  B0,  the  RHC'P  wave 
intersects  the  beam  line  two  times  as  indicated  by  temporal  growth  rates  at  two  separate 
regions  of  wave  number  k  for  a  given  axial  magnetic  field.  The  growth  rate  and  the  range  of 
unstable  values  of  k  increases  as  Ka  approaches  A'crit,3-  The  stability  boundaries  predict 
that  the  three-wave  instability  is  present  for  KCru,3  =  0.35  cm-1  <  K0  <  Kent, 4  = 
0.41  cm"1  (2.95  kG  <  B„  <  3.5  kG).  Numerical  results  for  It,  =  1  kA  are  in  excellent 
agreement  with  the  theory. 

Figure  7  plots  the  temporal  growth  rate  versus  k  in  the  orbit  unstable  region  for 
B0  =  2.15  kG  ( K0  =  0.26  cm-1).  The  region  of  instability  covers  both  positive  and 
negative  values  of  wave  number  k,  even  at  regions  where  the  beam  lines  are  far  from  the 
waveguide  mode.  Thus,  the  unstable  growth  rate  is  the  result  of  unstable  beam  orbits  in 
the  stellarator  and  guide  fields,  irrespective  of  the  electromagnetic  waves. 

As  the  beam  current  increases,  the  coupling  of  the  RHCP  and  LHCP  waves  becomes 
stronger  and  the  actual  stability  boundaries  deviate  from  the  analytic  expressions.  Figures 
8-10  are  plots  of  temporal  growth  rate  as  a  function  of  wave  number  k  for  Ib  =  10  kA. 
Figure  8  plots  the  temporal  growth  rates  in  region  I,  for  Ba  =  -1.0  kG  (Ka  =  -0.12  cm"1 ). 
B0  =  0,  Bo  =  1-0  kG  ( K0  =  0.12  cm-1),  and  Ba  =  1.3  kG  (Ka  =  0.156  cm-1). 

Figure  9  plots  temporal  growth  rates  versus  k  in  the  three-wave  unstable  region  II, 
for  B0  =  3.0  kG  ( K0  =  0.36  cm-1),  B0  —  3.25  kG  (K0  =  0.39  cm-1),  B0  =  3.5  kG  (K0  = 
0.41  cm"1)  and  Ba  =  3.85  kG  ( Ka  =  0.45  cm-1).  The  range  of  B„  for  the  three- wave 
unstable  region  II  at  zero  beam  current  is  2.95  kG  <  Ba  <  3.5  kG.  For  current  of 
Ip  =  10  kG,  the  RHCP  and  LHCP  waves  become  coupled  for  B0  =  2.95  -  3.25  kG. 
At  the  same  time,  the  three-wave  unstable  regime  is  extended  beyond  Ba  -  3.5  kG  to 
B0  =  3.85  kG.  The  RHCP  and  LHCP  waves  become  coupled  for  Ba  =  2.95  -  3.25  kG  and 
the  three-wave  unstable  regime  is  extended  into  a  region  where  Ka  >  Kct-ha  =  0.41  cm  -1 . 
where  the  analytic  (It,  =0)  expressions  predict  stability. 

Figure  10  plots  the  temporal  growth  rate  versus  k  in  the  orbit  unstable  region  for 
B0  =  2.15  kG  (K0  =  0.26  cm-1).  Here  the  region  of  instability  covers  essentially  all  values 
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of  wave  number  k. 

A  summary  of  peak  temporal  growth  rates  as  a  function  of  normalized  guide  field  I\0 
for  /&  =  10  kA  is  shown  in  Fig.  11  for  70  =  5,  y0  =  7.5  and  =  10.  The  group  of  curves 
on  the  left  belong  to  unstable  region  I  and  the  group  of  curves  on  the  right  belong  to 
unstable  region  II.  The  gap  separating  the  two  groups  of  curves  corresponds  to  the  orbit 
unstable  region.  Plots  of  the  group  velocity  as  a  function  of  Ka  are  shown  in  Fig.  12. 
Those  values  plotted  are  associated  with  the  growth  rates  of  Fig.  11.  The  group  velocity 
approaches  the  beam  velocity  in  three- wave  unstable  region  II. 

We  also  examine  the  effect  of  the  waveguide  radius  rg  on  the  three-wave  instability. 
As  the  guide  radius  is  increased,  the  waveguide  cut-off  frequency  unc  decreases.  For  the 
three-wave  unstable  region  I,  the  intersection  of  the  waveguide  mode  and  one  of  the  beam 
lines  in  Eq.  (28)  can  occur  for  negative  values  of  wave  number  k.  When  this  occurs,  the 
phaoe  velocity  is  negative  and  group  velocity  is  positive  but  reduced  in  value.  Figure  13 
plots  the  temporal  growth  rate  and  the  group  velocity  as  a  function  of  K„  for  rg  —  5  cm 
and  =  10  kA,  while  keeping  ail  the  other  parameters  the  same  as  in  Table  I. 

Comparison  of  Analytic  and  Numerical  Results 

The  analytical  expressions  for  the  temporal  growth  rates  for  the  two  regions  of  the 
three-wave  instability  are  given  in  Eq.  (42)  and  (44).  Figure  14  is  a  comparison  of  the 
maximum  temporal  growth  rate  as  a  function  of  the  normalized  guide  field  Ka  for  the 
values  obtained  from  the  numerically  solved  full  dispersion  relation  (solid  curves)  and 
from  the  analytical  expressions  (dashed  curves)  for  current  If,  =  1  kA.  The  agreement 
between  the  analytical  and  numerical  results  of  the  temporal  growth  rates  is  good  as 
long , aa  Au>2/c  >  T/c  (see  Eqs.  (39)-(44)).  On  the  dispersion  diagram,  2Au>;/c  is  the 
difference  in  frequency/c  of  the  beam  lines  (28)  that  could  go  unstable.  To  illustrate  the 
case  Au^/c  >  T/c,  we  take  Ka  =  —0.36  cm-1  (Ba  =  —3  kG).  The  numerical  result  of  the 
temporal  growth  rate  is  T/c  =  4.3  x  10-3  cm-1,  which  is  much  smaller  than  Au>2/c  =  0.25. 
The  analytical  result  for  the  temporal  growth  rate  is  also  T/c  =  4.3 x  10-3  cm-'.  As  A u>;;  c 
decreases  and  T/c  increases,  the  analytical  expressions  for  the  growth  rate  becomes  less 
accurate.  At  K0  =  U.12  cm-1  (B„  =  1  kG),  Awj/c  =  8.2  x  10-2,  which  is  comparable  to 
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h  rate,  V/c  =  2.50  x  10  2  cm  1 ,  is  8%  larger 


2  cm  1 . 


V.  CONCLUSIONS 


The  addition  of  stellarator  windings  to  an  axial  guide  field  was  proposed  as  a  method 
of  transporting  a  high  current  beam  in  a  curved  geometry  with  a  high  tolerance  to  energy 
mismatch.1  The  stabilty  properties  of  such  configurations  have  been  clarified  in  the  present 
study,  which  included  the  beam  centroid  motion,  the  electromagnetic  waveguide  modes, 
expressed  in  terms  of  right-hand  circularly  polarized  (RHCP)'  and  left-hand  circularly 
polarized  (LHCP)  waves,  and  the  induced  image  forces  on  the  conducting  boundaries,  in 
addition  to  the  external  magnetic  fields. 

We  find  five  operating  regimes:  i)  two  physically  distinct  three- wave  unstable  regimes, 
ii)  an  orbit  unstable  regime,  in  which  the  beam  centroid  is  unstable  independent  of  the 
electromagnetic  waveguide  modes  and  iii)  two  stable  regimes.  We  have  obtained  analytical 
expressions  for  the  boundaries  of  the  various  stability  regimes  in  parameter  space  and  have 
presented  algebraic  expressions  for  the  growth  rates  in  each  of  the  two  three- wave  unstable 
regimes.  These  analytical  results  are  valid  in  the  limit  of  low  beam  current,  when  the 
RHCP  and  LHCP  waves  decouple,  and  are  in  good  agreement  with  those  obtained  via 
numerical  solutions  of  the  full  dispersion  relation. 

The  simplified  dispersion  relation  for  uncoupled  RHCP  and  LHCP  waves,  given  in 
Eq.  (40),  has  been  shown  to  produce  results  that  are  in  close  agreement  with  those  of  the 
full  dispersion  relation.  The  simplified  dispersion  relation  is  not  valid,  however,  when  the 
coupling  between  the  RHCP  and  LHCP  waves  is  strong,  as  in  the  orbit  unstable  regime.  In 
three- wave  unstable  region  II,  the  coupling  is  strong  only  at  high  current.  The  growth  rates 
in  this  regime,  for  example,  from  the  full  (coupled)  dispersion  relation  are  about  20  —  30% 
larger  than  from  the  uncoupled  dispersion  relation  for  kw  =  0.5  cm~!  and  /  =  10  kA.  In 
three-wave  unstable  region  I  the  RHCP  and  LHCP  waves  are  essentially  uncoupled. 

These  results  suggest  that  the  three-wave  instability  can  be  avoided  by  appropriately 
choosing  the  various  parameters.  Results  show  that  as  the  beam  energy  increases,  the 
stability  conditions  become  more  restrictive.  In  such  cases,  it  may  be  necessary  to  decrease 
the  quadrupole  gradient,  Bqkq,  quadrupole  wave  number  kq  and/or  increase  the  value  of 
the  guide  field  B0  in  order  to  remain  in  the  stable  regime. 
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Table  I:  Parameters  Used  iu  Section  IV 


Quadrupole  gradient,  Bqkq 
Beam  energy  (ja) 

Beam  current, 

Drift-tube  radius,  rg 

Calculated  parameters 

TEXX  cut-off  frequency/c,  nu 

qnn 

Kqkq 

fio  =  Va/C 

For  rotating  quadrupole  wavelength.  A7  =  12.57 

Wave  number,  kq 
Kq 

Kcr-it, 2  {Bo  ~  1-32  kG) 

Kent, 3  [Bo  ~  2.95  kG) 

Kent ,4  (Bo  =  3.52  kG) 
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200  G/cm 
5 

1,10  kA 
3  cm 

0.614  cm'1 
0.256  cm~i 
0.024  cm~ 2 
0.9793 

0.5  cm'1 
0.048  cm'1 
0.154  cm"'1 
0.346  cm-1 
0.413  cm  - 1 
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'Plie  geometry  showing  the  mudurting  cylindrical  drift  tube,  the  rotating  quadrupnle 
Held  and  the  axial  guide  held. 


Wave  Number,  k  (cm  ') 

Fig.  2  The  dispersion  diagram  in  the  limit  of  zero  beam  current.  The  curves  are  the  LH 
and  RH  circularly  polarized  TEn  waveguide  modes.  The  dashed  (-  -  -)  and  solid  ( — ) 
straight  lines  are  the  beam  modes. 


Temporal  Growth  Rate/c,  f/c 


•ME 


Wave  Number,  k  (cm  1 ) 

Fig.  5  Plots  of  the  temporal  growth  rate  versus  k  in  three- wave  unstable  region  I.  for  current 
h  =  1  kA  and 

a)  B0  =  -1.0  k G  ( K0  =  -0.12  cm'1), 

b)  B0  =  0, 

c)  Bo  =  1.0  kG  ( K0  —  0.12  cm'1 )  and 

d)  B0  —  1.3  kG  ( K0  =  0.156  cm-1). 
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Temporal  Growth  Rate/c, 


2.95  kG 


0.020 


3.0  kG 


0.015 


0.010 


3.25  kG  2.95  kG 

/  3.0  kG  / 

3.5  kG  3.25  kG  \ 


0.005 


Wave  Number  k  (cm  l) 


Fig.  6  Plots  of  the  temporal  growth  rate  versus  k  in  the  three- wave  unstable  region  II.  lot 
current  /j,  =  1  kA  and 

a)  Ba  =  2.95  kG  ( K„  =  0.35  cm -1 ), 

b)  B0  =  3.0  kG  (K„  =  0.36  cm-1), 

c)  Ba  =  3.25  kG  (K0  =  0.39  cm-1)  and 

d)  S0  =  3.5  kG  ( K0  -  0.41  cm-1). 


Temporal  Growth  Rate/c,  T/c  (cm 


Wave  Number,  k  (cm  ') 


Fig.  7  A  plot  of  the  temporal  growth  rate  versus  k  in  the  orbit  unstable  region,  for  current 
h  =  1  kA  and  B0  =  2.15  kG  (K0  =  0.26  cm-1). 


Temporal  Growth  Rate/c,  17 c  (cm 


Fig.  9  Plots  of  the  temporal  growth  rate  versus  k  in  the  three- wave  unstable  region  II.  for 
current  /(,  =  10  kA  and 

a)  Ba  =  3.0  kG  (Ka  =  0.36  cm'1), 

b)  B0  =  3.25  kG  (K0  =  0.39  cm"1), 

c)  B0  —  3.5  kG  (K„  =  0.41  cm'1)  and 

d)  B0  =  3.85  kG  (Ka  =  0.45  cm~x). 
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Wave  Number,  k(cm  ') 

Fig.  10  A  plot  of  the  temporal  growth  rate  versus  k  in  the  orbit  unstable  region,  for  current 
lb  =  10  kA  and  B„  -  2.15  kG  ( Ka  =  0.26  cm'1). 


Temporal  Growth  Rate/c,  T/c  (cm  ') 


K0  (cm-1) 

Fig.  13  Plots  of  the  temporal  growth  rate  and  the  group  velocity  as  a  function  of  A'„  for 
rg  =  5  cm  and  current  =  10  kA,  while  keeping.  Other  parameters  are  the  same  as 
in  Table  I. 
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Group  Velocity/c, 


Temporal  Growth  R 


Fig.  14  Comparison  of  the  maximum  temporal  growth  rate  as  a  function  of  K„  for  the  values 
obtained  from  the  numerically  solved  full  dispersion  relation  (solid  curves)  and  from 
the  analytical  expressions  (dashed  curves)  for  the  parameters  in  Table  I  with  current 
ifc  =  1  kA. 
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Gaithersburg,  Maryland  20760 
Attn:  Dr.  Mark.  Wilson 

Naval  Postgraduate  School 
Physics  Department  (Code  61) 
Monterey,  CA  93940 
Attn:  Prof.  John  R.  Neighbours 
Prof.  Fred  Buskirk 
Prof.  Kai  Woehler 
Prof.  Xavier  Maruyama 

Naval  Surface  Warfare  Center 
White  Oak  Laboratory 
Code  R-41 

Silver  Spring,  Maryland  20903-5000 
Attn:  Mr.  W.  M.  Hinckley 
Dr.  M.  H.  Cha 
Dr.  H.  S.  Uhm 
Dr.  R.  Fiorito 
Dr.  K.  T.  Nguyen 
Dr.  R.  Stark 
Dr.  H.  C.  Chen 
Dr.  D.  Rule 
Dr.  Matt  Brown 
Mrs.  Carolyn  Fisher  (G42) 

Dr.  Eugene  E.  Nolting  (H23) 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 
Attn:  Dr.  C.  W.  Roberson 
Dr.  F.  Saalfeld 

Office  of  Naval  Research  (2  copies) 
Department  of  the  Navy 
Code  01231C 
Arlington,  VA  22217 


Princeton  University 
Plasma  Physics  Laboratory 
Princeton,  NJ  08540 
Attn:  Dr.  Francis  Perkins,  Jr. 

Pulse  Sciences,  Inc. 

600  McCormack  Street 
San  Leandro,  CA  94577 
Attn:  Dr.  Sidney  Putnam 

Pulse  Sciences,  Inc. 

2001  Vilshire  Boulevard 
Suite  600 

Santa  Monica,  CA  90403 
Attn:  Dr.  John  R.  Bayless 
Dr.  R.  Adler 

The  Rand  Corporation 
2100  M  Street,  NV 
Washington,  DC  20037 
Attn:  Dr.  Nikita  Wells 
Mr.  Simon  Kassel 

Sandia  National  Laboratory 
Albuquerque,  NM  87115 
Attn:  Dr. 
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David  Hasti/1272 

Collins  Clark 

John  Freeman/1241 

Charles  Frost 

George  Kamin/1274 

Gordon  T.  Leifeste 

Gerald  N.  Hays 

Michael  G.  Mazarakis/1272 

John  Wagner/1241 

Ron  Lipinski/1274 

James  Poukey 

Milton  J.  Clauser/1261 

Kenneth  R.  Prestwich/1240 

Kevin  O'Brien 

Isaac  R.  Shokair 
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Science  Applications  Inti.  Corp. 

5150  El  Camino  Road 
Los  Altos,  CA  94022 
Attn:  Dr.  R.  R.  Johnston 
Dr.  Leon  Feinstein 
Dr.  Douglas  Keeley 
Dr.  E.  Roland  Parkinson 

Science  Applications  Inti.  Corp. 

1710  Goodridge  Drive 
McLean,  VA  22102 
Attn:  Mr.  V.  Chadsey 
Dr.  A  Drobot 
Dr.  K.  Papadopoulos 
Dr.  William  W.  Rienstra 
Dr.  Alan  J.  Toepfer 
Dr.  Alfred  Mondelli 
Dr.  D.  Chernin 
Dr.  R.  Tsang 

Science  Research  Laboratory,  Inc. 

1600  Wilson  Boulevard 
Suite  1200 
Arlington,  VA  22209 
Attn:  Dr.  Joseph  Mangano 
Dr.  Daniel  Birx 

Commander 

Space  &  Naval  Varfare  Systems  Command 
PMV-145 

Washington,  DC  20363-5100 
Attn:  CAPT  J.  D.  Fontana 
LT  Fritchie 

SRI  International 
PS0-15 

Molecular  Physics  Laboratory 
333  Ravensvood  Avenue 
Menlo  Park,  CA  94025 
Attn:  Dr.  Donald  Eckstrom 

Dr.  Kenneth  R.  Stalder 

Strategic  Defense  Initiative  Org. 

SDI0/T/DE0 

The  Pentagon 

Washington,  DC  20009-7100 
Attn:  Lt  Col  R.  L.  Gullickson 
Dr.  D.  Dus  ton 

Titan/Spectron,  Inc. 

P.  Q.  Box  4399 
Albuquerque,  NM  87196 
Attn:  Dr.  R.  Bruce  Miller 
Dr.  John  Smith 


Titan  Systems,  Inc. 

2685  Marine  Way 
Suite  1408 

Mountain  View,  CA  94043 
Attn:  Dr.  Kenneth  U.  Billman 

Titan  Systems,  Inc. 

9191  Tovne  Centre  Dr. -Suite  500 
San  Diego,  CA  92122 
Attn:  Dr.  R.  M.  Dove 

University  of  California 
Physics  Department 
Irvine,  CA  92664 
Attn:  Dr.  Gregory  Benford 
Dr.  Norman  Rostoker 

University  of  California 

San  Diego,  CA  92110 

Attn:  Dr.  Marshall  N.  Rosenbluth 

University  of  Maryland 
Physics  Department 
College  Park,  MD  20742 
Attn:  Dr.  Y.  C.  Lee 
Dr.  C.  Grebogi 
Dr.  W.  Destler 
Dr.  C.  Striff-ler 

University  of  Michigan 
Dept,  of  Nuclear  Engineering 
Ann  Arbor,  MI  48109 
Attn:  Prof.  Terry  Kammash 
Prof.  R.  Gilgenbach 

Director  of  Research 
U.S.  Naval  Academy  - 
Annapolis,  MD  21402  (2  copies) 


Do  NOT  make  labels 
for  these  tivo--b e  1  oui: 
Records - ( 1  copy  ) 

Naval  Research  Laboratory 
Washington,  DC  20375-5000 
Code  2630 
Timothy  Calderu/ood 


